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Abstract. We present a variational principle governing the quasistatic evolu- 
tion of a linearized elastoplastic material. In case of linear hardening, the novel 
characterization allows to recover and partly extend some known results and 
proves itself to be especially well-suited for discussing general approximation 
and convergence issues. In particular, the variational principle is exploited 
in order to prove in a novel setting the convergence of time and space-time 
discretizations as well as to provide some possible a posteriori error control. 



1. Introduction 

The primal initial-boundary value problem of clastoplasticity consists in deter- 
mining the generalized deformation state of a material subject to external me- 
chanical actions. In particular, starting from some initial state and for a given load 
and traction, one shall determine the displacement u of the body from the reference 
configuration, the inelastic (plastic) part p of its strain, and, possibly, a vector of in- 
ternal hardening variables ^. In the small deformation regime and within the frame 
of associative elastoplasticity, the problem is classically formulated in a variational 
form as that of finding the absolutely continuous trajectory t G [0,T] i-^ y{t) G Y 
{Y is a Banach space) such that 

diiy')+Ay3e a.e. in (0,T), 2/(0) - yo, (1-1) 

where y = {u,p,(,) stands for the vector of unknown fields, A : Y ^ Y* (dual) 
is linear, continuous, and symmetric, and ip '■ Y [0, oo] is the positively 1- 
homogeneous and convex dissipation potential {d is the classical subdifferential in 
the sense of Convex Analysis, see below). Moreover, £ : [0,T] — > F* is a given 
and suitably smooth generalized load (possibly including surface tractions) and 
yo represents the initial state. The reader is referred to Section [5] for some brief 
mechanical motivation as well as to the classical monographs by DuvAUT & LiONS 
[DL76| . Han & Reddy |HR99j . Lemaitre & Chaboche |LC90j . and Simo & 
Hughes |SH98| for a comprehensive collection of results. 

The aim of this paper is that of investigating a global-in-time variational formula- 
tion of problem (jl.ip . In particular, we shall introduce the functional 
T : W^^^{0,T;Y) [0,oo] on trajectories as 

Hy) = r W) + r{i- Ay) - (£ - Ay, y') ) , 
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where ip* stands for the conjugate tp*{w) = sup^,gy ((w, w) — ^p{v)) of ip and (•,•) 
denotes the duahty pairing between Y* and Y. The starting point of this analysis 
relies on the fact that solutions of (jl.ip and minimizers of fulfilling the initial 
condition coincide, namely (see Theorem 13. 1|1 

y solves dLlI) iff T{y) =mmT = and y{0) = y„. (1.2) 

This variational characterization has a clear mechanical interpretation. Indeed, 
since '0 is positively 1-homogeneous, its conjugate ip* turns out to be the indi- 
cator function of the convex set dip{0). Hence, minJF is actually a constrained 
minimization problem and we have 

e-Aye dipiO) a.e. in (0,r), 

Hy)=0 iff { ^i^T,y{T))+ r i:{y')^^{Q,ym- r {t'^y), ^^'^^ 







where we have used the notation (t, y) ^ (p{t,y) = i {Ay,y) — {£{t),y). The first 
relation above expresses the so-called local stability |Mie05j of the trajectory whereas 
the second is nothing but the energy balance at time T. More precisely, ip{t,y) 
denotes the complementary energy at time t for the state y, /J" "0(2/') represents the 
dissipation of the system on [0, T], and — jj" y) is the work of external actions 
on [0,r]. Hence, minimizing T consists in selecting the (only) stable trajectory 
which conserves the energy. In this regards, the reader is referred to the pioneering 
papers by MoREAU [Mor68( iMorTOl IMor71| . 

The interest of variational characterization (|1.2|) of the differential problem (jl.ip 
relies on the possibility of exploiting the general tools from the Calculus of Vari- 
ations. Some care is however required. Indeed, although J- is convex and lower 
semicontinuous with respect to the weak topology of iy^'^(0, T; F), the functional 
generally fails to be coercive. Moreover, one is not just asked to minimize T but 
also to prove that the minimum is 0. This considerations suggest that the Direct 
Method is hardly applicable in order to get solutions to (|l.ip via the characteriza- 
tion in ([Ll]). 

The first issue of this paper is instead that of showing that the variational prin- 
ciple in ()1.2p is particularly well-suited for discussing general approximation issues. 
Since solutions and minimizers coincide, a natural tool in order to frame an abstract 
approach to limiting procedures within (jl.ip is that of considering the corresponding 
minimum problems via F-convergence |GF75j . As the value of the functional is di- 
rectly quantified to be on the minimizers, what is actually needed here for passing 
to limits are so-called F-lim inf inequalities only and the latter are generally easily 
available. We shall specifically focus on the case of linear hardening elastoplasticity 
and apply the above-mentioned perspective in order to recover in a unified and 
more transparent frame and partly generalize some convergence results for confor- 
mal finite elements (Thm. 15. 3p . time-discretizations (Thm. 16. 5p . and fully-discrete 
space-time approximations (Thm. 17. ip . In particular, for time-discretization we 
develop a discrete version of the variational principle ()1.2p in the same spirit of 
the theory of variational integrators [MWOlj (see Subsection 16. ip . This connec- 
tion entails also some generalized view at the classical discrete time-schemes (see 
Subsection 16. 5p . 

A second novel point of the present variational approach consists in the possibility 
of exploiting J- in order to estimate a posteriori some approximation error. By 
letting T{y) = 0, we will check that (Cor. 14. 5p 

Yia.^\{A{y~v),y-v) <T{v) V« e W^^'^O, T; F), i;(0) = j/o- 
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If A shows some coercivity (which is precisely the case of hnearized hardening, 
see Subsection 12 ■4p . and v is the outcome of some approximation procedure, the 
estimate above may serve as the basis for some a posteriori estimation procedure, 
possibly headed to adaptivity (see Subsection l6.7p . Let us stress that the latter and 
(II. 3p entail that the distance of a (stable) trajectory from the solution to (II. ip can 
be uniformly estimated by means of its energy production along the path. 

The variational characterization in (|1.2p stems from a reintcrpretation in the 
present elastoplastic frame of the celebrated principle by Brezis & Elekand 
jBE76al [i3E76b| for gradient flows of convex functional. Since its introduction, the 
latter principle has continuously attracted attention. In particular, it has been ex- 
ploited in the direction of proving existence |Rio76a[ IRio76b[ IRio78i IRio79i IAuc93[ 
IRou00|[GT04| (note that the above-mentioned obstructions to the application of the 
Direct Methods again appear) and the description of long-time dynamics |Lem96| . 
Moreover, the Brezis- Ekeland approach has been adapted to the case of second 
order |Mab01i IMab03j and doubly nonlinear equations |Ste06aj as well. 

One has to mention that, of course, (|1.2p is not the only possible global-in-timc 
variational characterization of (|l.ip . Besides minimizing the space-time norm of 
the residual (which might be little interesting since the order of the problem is dou- 
bled) , one has at least to mention Visintin jVisOlj , where generalized solutions are 
obtained as minimal elements of a certain partial-order relation on the trajectories, 
and the recent contribution by Mielke & Ortiz jMO07| where the functional 

y ^ e-^/^m y(T)) + £ e-*/- (^V(y') + -^^{t, y)) (1.4) 

is minimized among trajectories with ?;(0) = yo. Under extra-smoothness conditions 
on Tp (not fulfilled in the current frame) , the Euler-Lagrange equations of the latter 
functional are 

-eDV(y')y" + D^Piy') + Ay = £ 
y(0) = yo, D^iy'iT)) + Ay{T)=e{T). 

In particular, minimizing the functional in (|1.4p consists in performing a suitable 
elliptic (in time) regularization of the problem. In the specific case of ip positively 
1-homogeneous, the limit e — > can be carried out and the minimizers of the 
functional in p.4p are proved to converge to the solution of p.ip . The latter 
approach is quite different form that of (|1.2p . On the one hand, it is much more 
general as it naturally applies to the non-smooth case as well (no derivatives of tp 
and (/) are involved). On the other hand, it relies on an intermediate and somehow 
unphysical (not causal) e-regularized problem and (besides existence) it seems not 
directly suited for recovering the full extent of our approximation results for the 
specific case of problem (jl.ip . 

2. Mechanical model 

Let us provide the reader with a brief introduction to the mechanical setting 
under consideration. Our aim is just that of recalling some essential features of 
the models and well as their variational formulation. In particular, we restrain 
from reporting here an extensive discussion on associative elastoplasticity as the 
latter can be easily recovered from the many contributions on the subject. The 
reader is particularly referred to the mentioned monographs for some comprehensive 
presentation. 
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2.1. Preliminaries. We will denote by M^y^ the space of symmetric 3x3 tensors 
endowed with the natural scalar product a : b := tr(a6) — aijbij (summation 
convention). The space Kgym is orthogonally decomposed as M^y^-, = R^^^ ©IRI2, 
where M I2 is the subspace spanned by the identity 2-tensor I2 and IR'^g^ is the 
subspace of deviatoric symmetric 3x3 tensors. In particular, for all a G Kgymi we 
have that a = a^cv + tr(a)l2/3. 

We shall assume the reference configuration 51 to be a non-empty, bounded, and 
connected open set in R'^ with Lipschitz continuous boundary. The space dimension 
3 plays essentially no role throughout the analysis and we would be in the position of 
reformulating our results in with no particular intricacy. Our unknown variables 
are the displacement of the body u S M'^, the plastic strain p € Rjov' ^^"^ ^ vector 
of internal variables ^ G M™ (to G N) which will describe the hardening of the 
material. We will denote by e{u) the standard symmetric gradient. 

2.2. Constitutive relation. Moving within the small-strain regime, we additively 
decompose the linearized deformation e{u) into the elastic strain e and the inelastic 
(or plastic) strain p as 

e(u) = e + p. 



five definite linear map C : R^y^, I^symi shall assume that the orthogonal 
subspaces R'^^^ and RI2 are invariant under C. This amounts to say that indeed 



Let C be the elasticity tensor. By regarding the latter as a symmetric posi- 

3x3 _^ K3x3 
sym sym ' 

ID) 1 _ W 

dcv 

Ca ~ Cdcvfldcv + Ktr(a)l2, 

for a given Cdcv : I^dov ~^ I^dev ^^'^ ^ constant k, and all a G R^y^. The case 
of isotropic materials is given by Cdev = 20(14 — I2 <E) I2/3) and G and k are 
respectively the shear and the bulk moduli. The latter decomposition is not ex- 
ploited in our analysis but it is clearly suggested by the mechanical application. 
Moreover, we shall introduce two linear symmetric positive semi-definite harden- 
ing moduli Hp : R^^^^ R^^^^ and : R"" ^ R™ (to be identified with a 
fourth order tensor and a matrix, respectively) and define the Helmholtz free en- 
ergy W : R3^„3 X M^j^^ X R" ^ [0, 00) of the material as 

W{e{u),p, ■■= l{e{u) - p) : C{e{u) - p) + ^p : Mpp + ■ 

The generalized stresses (cr, ry) arc conjugate to the above-defined generalized 
strains (e,^) via the energy W. In particular, the material is classically assumed 
to show elastic response, 

dW 

a= — = Ce^C{e{u)-p). (2.1) 

and the thermodynamic force 77 driving the evolution of the internal variables ^ is 
defined as 

dW 

^/ = -^ = -%e. (2.2) 

Moreover, moving within the frame of associative elastoplasticity, we assume the 
existence of a function R : R^^^ x R™ — + [0,oo] convex, positively 1-homogeneous, 
and lower semicontinuous such that 

^Rip,0^(;'-^^'). (2.3) 

In particular, R is asked to be the support function of a convex set C* G Rdov ^ 
i.e. R{p) = supggp, q : p. We will indicate with R* its conjugate, namely the 
indicator function of C* given by R*{q) = if q G C* and R*{q) ~ 00 otherwise. 
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Moreover, we let C be the domain of R, namely C = D(R) ^ {{p, ^) e M.^^^ x R™ : 
R{p,^) < oo}. 

Finally, the above material relations (|2.1|l - (|2.3p can be condensed as the following 
constitutive material law 

a«(«,.(<-«''>n">). 

which in turn can be rephrased in the form of (jl.ip by letting 

y = y = M'li^ev X K'" , i' = R, 

= ((C + H>,%e), i={Ceiu),0). (2.5) 

Let us close this subsection by explicitly mentioning three classical linear hardening 
models |HR991 Ex. 4.8, p.88] 

Linear kinetic hardening: choose Hp = hpl4 where hp > 0, and = 0. 

In this case the internal variable ^ is not evolving and shall be removed 

from the set of unknowns. 
Linear isotropic hardening: choose Hp = Q, m ~ 1, and H^ = ft,^ > 0. 

Moreover, let D{R) = {(p, G K^^^ x M : |p| < £,}. 
Linear combined kinetic-isotropic hardening: let Hp = /ipl4, m = 1, 

and = where hp, h^ > 0. Moreover, let D{R) = {{p,^) e R^^^^ x M : 

ipi<a- 

It is beyond the purpose if this introduction to discuss and justify the above- 
mentioned material models. The reader shall check the cited references for com- 
ments on their relevance within applications and some mechanical motivation. 



2.3. Variational formulation of the quasistatic evolution. Let us now move 
to the consideration of the full equilibrium problem. To this aim, we assume that the 
boundary dQ is partitioned in two disjoint open sets Ftr and Foir with dVtr = 9FDir 
(in dQ). We ask Foir to be such that there exists a positive constant CKorn depending 
on Foir and Q such that the Korn inequality 

CKorn||lt||Hi(0;R3) < || « || i2(rj,„.R3) + II ^2 (q.rS x 3 ) , (2.6) 

holds true for all u £ iJ^(ri;R'^). It would indeed suffice to impose Foir to have a 
positive surface measure (see, e.g., |DL761 Thm. 3.1, p. 110]). 

For the sake of simplicity, we will prescribe homogeneous Dirichlet boundary 
conditions on Foir (our analysis extends with little notational intricacy to the case 
of non- homogeneous Dirichlet boundary conditions as well). On Ftr some time- 
dependent traction will be prescribed instead. 

As for the full quasistatic evolution of the material we shall couple the constitu- 
tive relation ()2.4p with the equilibrium equation 

divCT + / = in Q. (2.7) 

Here, we assume to be given the body force / : [0, T] iy^(il;]R'^) and a surface 
traction g : [0,T] ^ L2(Ftr;K3). 
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Then, one can rephrase the problem into the form of (|l.ip by choosing 



(2.8) 



such that u = Q on Fdi 



(2.9) 



{A{u,p,£,),{v,q,z)) = I ({e{u)~p):£{e{v)-q)+p:mpq + e -^^z 




(2.11) 



(2.10) 



and defining the total load i : [0, T] y* as 



(^(t),(u,p,e)) = / f-u + [ g-udn^ ViiG i/^(17;R3), t e [0,T], 



where Ti^ is the 2-dimensional HausdorfF measure. 

2.4. The coercivity of A. Let us close this introductory discussion by explicitly 
commenting on the coercivity of the bilinear form induced by A. We shall recall 
some sufficient conditions on Hp, H^, and R in such a way that there exists a 
constant a > such that 



where | • | is the norm in Y. This issue is fairly classical jHRQQ) Sec. 7.3, p. 167] 
and we discuss it here for the sake of completeness only. 

Of course (|2.12p holds (and even for all y € Y) whenever Hp and Hj are positive 
definite (this is the case of the above-mentioned linear combined kinematic-isotropic 
hardening). 

As we have already observed, in case H^ = 0, the problem naturally reduces to 
the pair {u,p) only. Up to this reduction, (|2.12p holds (again for all y £ Y) when 
Hp is positive definite. This is exactly the case of linear kinematic hardening. 
On the other hand, in case Hp = 0, the plastic strain will still evolve and one 



has ([21^ if D{R) is bounded in the p-direction for all ^, namely if |HR99[ (7.51)] 
D{R) C {{p,0 e KdcvXlI^" : /^bP < for some constant /3 > 0}, (2.13) 



which is clearly the case for linear isotropic hardening. 

Some generalization of the latter condition could in principle be considered for 
the case when Hp and Hj are only semi-definite. In particular, (|2.12p holds if one 
assumes (|2.13p and 



Let us mention that the most critical case in the class of (|2.4p is Hp = 0, 
Hj = where actually no hardening takes place. This is the situation perfect 
plasticity for which the Sobolev space framework above is not appropriate and one 
would consider the space BD{Vl) of functions of bounded deformations instead 
[DDM06] . We shall make clear that, even if our variational characterization covers 
the case of perfect plasticity, the subsequent approximation results apply to the 
linear hardening situation only. 



{Ay,y)>a\y\'' ^y £ D{i,) 



(2.12) 



C^O and e'^.%e = ^ i?(p,0 



3x3 
■dcv ■ 
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3. Characterization 

3.1. General assumptions. Let us start by recalling notation and enlisting the 
basic assumptions for the following. First of all, we will ask that 

Y is a separable and reflexive Banach space. (3.1) 

We will use the symbols | • | for the norm of Y and (•, •) for the duality pairing 
between Y* (dual) and Y. The norm in Y* will be denoted by | • |* instead. 

We introduce the functional 

tp : Y [0, oo] proper, convex, lower semicontinuous, 

and positively 1-homogeneous. (3.2) 

Equivalently, ip is required to be the support function of a convex and closed set 
C* C Y* containing 0, namely 

^(y)=sup{(2/*,y) : y*eC*}. 

We shall define C = D{ip). Hence, the conjugate ip* : Y* ^ [0, oo], which is 
classically defined as ip*{y*) = sup,ygy((2/*, y) — ^{y))^ is the indicator function of 
the convex set C*, namely 4>*{y*) = if y* e C* and '4>*{y*) = oo otherwise. Let us 
remark that ip fulfills the triangle inequality ipi'^) ^ V'(^) +'0(c) whenever a = b + c. 

We shall use the symbol d in order to denote the usual subdifferential in the 
sense of Convex Analysis, namely 

y* G dipiy) iff y & D{ijj) and (y*, w ~ y) < ip{w) — ip{y) Vw g Y. 

Similarly, we define 

yedip*{y*) iff y* D{i]f) and {w* ~ y\y) < i:* {w*) - t^* {y*) ^w* (^Y* 

iff y* e C* and {w* ~ y\y) < Vw* e C* . 
Finally, we recall Fenchel's inequality 

^P{y) + r{v*)>{v\y) yyeY,y*eY*, 

and remark that equality holds iff y* G dip{y) (or, equivalently, y G d*tp*{y*))- 
As for the operator A we require 

A : Y ^ Y* linear, continuous, and symmetric, (3.3) 
and define the function 

y Hy) = 2 i^^y^y) > 

so that A = D<j>. Moreover, we will ask to be coercive on C = D{il'), namely we 
assume that there exists a positive constant a such that 

Hy)>^\y\' ^y^c. (3.4) 

As we have already commented in Subsection 12.41 the latter coercivity is fulfilled 
in the situation of elastoplastic evolution with linear kinematic, isotropic, or com- 
bined kinematic-isotropic hardening and will turn out to be sufficient for both the 
forthcoming characterization results. 

On the other hand, the following uniqueness-type results will be checked under 
some stronger coercivity frame and we will ask for 

0(y)>f|yP VyeC-C. (3.5) 

Clearly, condition (|3.5p is fulfilled when (f> happens to be coercive on the whole space 
Y. The latter applies in particular to the case of linear kinematic and combined 
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kinematic-isotropic hardening elastoplasticity. In this case, (f> defines an equivalent 
(squared) norm in Y. 

We shaU make use of the following notation 

x(y) = Hy) + Vy e r 

Indeed the latter choice is just motivated by simplicity and could be replaced as 
well by any other x ■ Y s- [0, 00) such that xiv) = iff y = and y i— > xiv) ~ 4'{y) 
is lower semicontinuous. 

Finally, we shall fix data such that 

ieL°°{0,T;Y*), yoGC. (3.6) 

The restriction on the choice of the initial datum in C is motivated by the coercivity 
assumption on in (|3.4p . On the other hand, we shall explicitly mention that the 
usual choice for yo hi elastoplasticity is yo = 0. 

In the forthcoming of the paper the above assumptions ()3.ip - (|3.4p and ()3.6p will 
be tacitly assumed (unless explicitly stated). It should be however clear that the 
above choice is motivated by the sake of simplicity. Indeed, most of the following 
results still hold under suitably weaker assumptions, as we shall comment. 

3.2. The functional. Let the Lagrangian i : (0, T) x y x y — > [0,oo] be defined 
as 

L{t, y,p) = ^{p) + rm - Ay) - (e{t) - Ay,p} 

for a.e. t € (0, T), Vy, peY, (3.7) 
and the functional F : W^'^{0,T;Y) -> [0,oo] as 

F(y) = / L{t, y{t), y'{t)) dt + x(y(0) - yo). (3.8) 

JO 

Now, by simply using the chain rule, we obtain that 

F{y) = ^ (V(y') + - Ay) - {I, y') ) + 0(y(r)) - 0(y(O)) + x(2/(0) - yo). 
A first remark is that, by exploiting the particular form of 
F(y) = (^^{y')+r{t-Ay)~{i,y')) 

+ 0(y(r))+0(yo)-(Ay(O),yo) + |y(O)-yop. (3.9) 
In particular, F is clearly convex. 

3.3. The characterization. 

Theorem 3.1 (Variational principle), y solves iff F{y) = = mini^. 

Proof. Owing to Fcnchcl's inequality we have that 

L{t, y,p)^Q ifr t{f) - Ay e 9V(p) 

and, clearly, x(2/(0) — yo) = iff y(0) = yo. Hence, all solutions y of (jl.ip are such 
that F{y) = and vice-versa. □ 

Let us remark that the latter variational characterization result holds in much 
grater generality. The proof made no use of the separability and reflexivity of Y 
nor of the linearity of A (besides its being single- valued and such that 1 1— > Ay{t) is 
measurable). Moreover, the positive 1-homogeneity of is unessential |Ste07| . In 
particular, the variational approach of Theorem 13. II can be directly extended to a 



A VARIATIONAL PRINCIPLE FOR HARDENING ELASTOPLASTICITY 



9 



variety of different dissipative systems possibly including viscous evolution as well. 
We shall address this perspective in a forthcoming contribution. 

We have already observed that F is convex. Moreover, F is lower semicontinuous 
with respect to the weak topology of VF^''^(0, T; Y) since all weakly convergent se- 
quences in W^'^{0, T; Y) are pointwisc weakly convergent as well. Hence, one could 
be tempted to use the Direct Method in order to get the existence of minimizers, 
i.e. solutions to equation p.ip . As we commented in the Introduction, this seems 
to be no trivial task. 

First of all, the functional F need not be coercive with respect to the weak 
topology of W^'^{0,T;Y). Indeed, the functional ip may degenerate and hence 
not control the norm of its argument. Moreover, even in the case when ip is non- 
degenerate, the homogeneity assumption just entails that the sublevels of F are 
bounded in W^'^{0,T;Y) and no weak compactness follows. 

Secondly, even assuming cocrcivity in the weak topology of W^'^{0,T;Y), one 
would still need to prove that the minimum is attained. This very obstruction to 
the use of the Direct Method occurs for the Brezis & Ekeland principle for gradient 
flows }BE76b| Rem. 1] and for its doubly nonlinear generalization in |Ste06aj as 
well. 



3.4. The variational principle for hardening elastoplasticity. By referring 
to the notations of Section [2 let us now present the actual form of the functional 
F for the case of the constitutive relation for linearized elastoplastic materials with 
linear hardening (see (|2.5p ). In this case the functional reads 



i(p(0) -po) : (C + MpMO) - po) + - Co)^ • He(C(0) - Co) 

|(p(0),C(0))-(po,Co)| 



2 



for some given initial datum (po,Co) e ^dcv e{u) e L°°(0, T; M^^^). 

In the situation of the quasistatic evolution, for some given initial datum 
(mo,Po,Co) e Y, a load / e L°°{0,T; {L^{n;R^))), and a traction 
g e L°°(0,T; {L'^(Ttr;M.^))), the functional reads (see (jTHI - lj?!!]) ) 

F{u,p,0 = ^ (i?(p,O + i?*(C(e(u)-p)-Hpp,-%0) 
[ f-u- [ [ g-um'' 

Jo JQ JTtr 

T 



Jn 



{e{u) - p) : C{e{u) - p) + p : H^p H- • H^c) 



i / (£(u(0) - ua) - (p(0) - po)) : C{e{u{Q) - u^) - (p(0) - po)) 

1 
2 

i / |(ii(0),p(0),e(0))-(7.o,Po,eo)P 



(p(0) - po) : Hp(p(0) - po) + {m - Co)^ • ~ Co) 
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for all points {u,p,^) G Y such that 



(£(u) - p) : Ce{v) = f -v- g-vdW 

\/v e H'^{n;M.^) with v = on Toir, a.c. in (0,r), 
and F(u,p,^) ~ oo otherwise. 

4. Properties of the minimizers 

For the sake of illustrating the variational principle of Theorem 13.11 we shall 
collect here some properties of the trajectories belonging to the domain of the 
functional F and, in particular, of the minimizers. 

4.1. Trajectories are in C. 

Lemma 4.1. Let F{y) < oo. Then y{t) e C for all t e [0,T]. 

Proof. Since F{y) < oo we have that y' G C almost everywhere in (0, T). Hence, 
for all t £ [0, T], we have that y' e C by Jensen's inequality. On the other hand 
yo £ C and y{t) = yo + Jq y' ■ The assertion follows by recalling that C is a cone. □ 

4.2. Stability at regular points of £. Assume £ : [0, T] Y* is given and let 
the set of stable states S{t) C Y ioi t E [0,T] be defined as 

S{t) ={yeY : 4>{y) - y) < (f>{w) ~ {£{t),w) + ^{w - y) Vu; G Y}. 

One can easily prove that 

y e S{t) iff e{t) -Aye 90(0) = C*. 

Lemma 4.2 (Stability of the minimizers). Let £ be either left- or right-weakly 
continuous at some point t £ [0, T] and F{y) < oo. Then y{t) £ S{t). 

Proof Since F{y) < oo, we have that e{t) - Ay{t) e= C* for aU t £ (0,T) \ iV, 
where \N\ = 0. Choose a sequence tk £ (0, T)\N such that tk t (from the left or 
from the right) and £{tk) e{t) weakly in Y*. Hence, £{tk)-Ay{tk) -> e{t)-Ay{t) 
weakly in Y* and e{t) - Ay{t) £ C* . □ 

In particular, if i happens to be right-continuous at 0, the functional F will not 
attain the minimum value unless the initial datum yo is stable, namely £ >S'(0). 

4.3. Equivalent formulations. Letting now I £ VF^'^(0, T; F*), problem (jl.ip 
admits some alternative equivalent formulations |Mie05[ Sec. 2.1]. We explicitly 
mention that y £ W^'^{Q,T;Y) is said to be an energetic solution if it solves the 
energetic formulation jMT04j of p.ip , namely 

y{t)£S{t) Vie[0,T], (4.1) 
q^{yit))-m,yit))+ [ V^y') = 0(2/(0)) - (^(0), y(0)) - / (£',y) 



Vie[0,T], (4.2) 
2/(0) = 2/0. (4.3) 

MiELKE & Theil [MT04j proved that the latter is equivalent to ()l.ip and hence, 
owing to the characterization of Theorem 13. 11 to F{y) = = miiiF. For the aim 
of pointing out some features of our variational approach, we shall present here a 
direct proof of this fact. 
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Lemma 4.3 (Equivalence with the energetic formulation). Let £ S W^'^{Q,T;Y*). 
Then, F{y) = = minF iff y fulfills t^TDf - ^ . 



Proof. Owing to Lemma we readily have that the stability condition (|4.ip holds 
iff ip*{£ — Ay) ~ almost everywhere. 

Let y be such that F{y) = 0. Then gU) and (gS]) hold and L(t, y{t),y'{t)) = 
for a.e. t e (0, T). In particular, for all t G [0, T], 



= Lis, yis), y'is)) ds = {^P{y') - {£, y') ) + 



so that the energy equality (|4.2p holds for all t G [0, T]. 

On the contrary, let y G ^^^'^(O, T; Y) fulfiU gl])-g31). Then x(y(0) - ?;o) = 
and ?/'*(^ — Ay) = almost everywhere (see above). Hence F{y) = follows from 
the energy equality (|4.2p at time T and an integration by parts. □ 



Let us mention that the latter lemma proves in particular that the energy equality 
(|4.2p could be equivalently enforced at the final time T only. Moreover, it proves 
that, as already commented in the Introduction, all stable trajectories t i— > y(t) 
(i.e. trajectories such that y{t) G S{t) for all t G [0,T]) are such that the following 
energy inequality holds 

Hyit))-{m,y{t))+ f 'W)><p{ym~m,ym~ f\e\y) vtG[o,T]. 

Jo Jo 

Hence, we have provided a proof to [Mie05|, Prop. 5.7] (in a somehow simpler 
situation though). 

Before closing this subsection, let us explicitly remark that the above inferred 
equivalence between formulations has been obtained for the absolutely continuous 
case £ G W^'^{0,T]Y*) only, whereas the characterization of Theorem 13.11 holds 
more generally for bounded £. 



4.4. The functional controls the uniform distance: uniqueness. So far, we 
have simply reformulated known results in a variational fashion. Here, we present 
some novel results instead. 



Lemma 4.4 (Uniform distance control via F). We have that 

77(1 - ri) rnax (j){u{t) - v{t)) < r]F{u) + (1 - r))F{v) 

Vm, w G W^^'^(0,T;y), ry G [0,1]. (4.4) 



Proof. The statement follows from the quadratic character of (j). Fix t G [0, T] and 
define F* : W^^^[0,t-Y) [0,oo] as 

F\y) = / L{s,y{s),y'{s))ds + x{y{^)-yo) 
Jo 

t 



iiy')+r{^-Ay)-{£,y')^ 
^{y{t)) + 0(j/o) - {Ay{Q),yo) + |?/(0) - yo? ■ 
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Then, clearly y i-+ G*(y) = F*{y) — 4'{y{t)) is convex. Hence, letting w = t]u + {1 — 
r])v we have that 

< F*{w) 

< ii{G\u) + cf>{u{t))) + (1 - v){G\v) + ^iv{t))) - ry(l - rj^uit) - v{t)), 
whence the assertion follows. □ 

The latter lemma exploits the quadratic character of (/> only. In particular, no 
coercivity for cj) is assumed. It should however be clear that its application (as the 
title of the lemma indeed suggests) will always be referred to the situation where 
the stronger (|3.5p is required, namely when the left-hand side of (|4.4p controls 

77(1 — 77)— max \u — wP. 

^ ^ 2 [0,T] ' ' 

We now present two immediate corollaries of Lemma 14.41 
Corollary 4.5 (Uniform distance from the minimizer). Let F{y) = 0. Then 

max (l)(y(t) - v(t)) < F{v). 

te[o,T] ^ ' 

This corollary encodes an interesting novel feature of our variational approach, 
for it provides possible a posteriori error estimator to be used within approximation 
procedures. It is interesting to remark that the uniform distance of any stable 
trajectory from the minimizer is controlled by means of its energy production along 
the path only. Again, although Corollarv l4. 51 holds under no coercivity assumptions 
of let us mention that its application will be restricted to the frame of (j3.5p . 
Finally, we have uniqueness of the minimizers of F attaining the value 0. 

Corollary 4.6 (Uniqueness). Assume p.5|) . Then, there exists at most one tra- 
jectory y such that F(y) = 0. 

4.5. Lipschitz bound. Form this point on and throughout the remainder of the 
paper we shall tacitly assume 

£e iyi.-(o,r;r*), yoeS{Q). (4.5) 

As already commented after Lemma 14.21 the above restriction on the initial 
datum is mandatory whenever i admits a weak-right-limit in 0. 

As for £, the extra Lipschitz continuity assumption is motivated by the rate- 
independence of the problem (every absolutely continuous datum can be time- 
rescalcd to a Lipschitz continuous datum) and the following well-known result. 

Lemma 4.7 (Lipschitz bound). Assume and let £ € W^°°{0,T;Y*), and 

F{y) = 0. Then 

II2/'IIl-(o,T;F) < ^II^'I1l-(o,T;Y.) a.e. m {0,T). (4.6) 

The proof of the lemma is exactly the classical one jMT04| Thm. 7.5], but 
formulated by means of our variational arguments. We provide it for the sake of 
completeness. 

Proof. Let 0<s<i<rbe fixed. Since L{y, y') = almost everywhere we have 
that 



(V'(y') + (^',y)) + 



0. 
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On the other hand, owing to the strong monotonicity of A and the fact that y{s) G 
S{s) (see Lemma one obtains 

Hvit) - y{s)) < Hvit)) - m,y{t)) + ^{y{t) - y{s)) - ^{y{s)) - {({s), y{s)) 
= (</.(y) - {£, u)) \ [ (/(r), y{t)) dr + i>{y{t) - y{s)). 

« J s 

By taking the sum of these two relations and rccalhng that, by Jensen, 
^P{y{t) - y{s)) = ^ ( I y'\ < f V(y') 



we get that 

<l>{y{t) - y{s)) < I {e'{r),y{t)-y{r)) dr. 



Finally, an application of some extended Gronwall lemma (see |Mie05|, Thm. 3.4]) 
entails that 

f l2/(t) - y{s)\' < l\\nL^io,T;Y')\y{t) ~ y{s)\ (t - .) 

and the assertion follows. 

□ 



5. Space and data approximation 

We now apply the characterization results of Theorem 13.11 to the approximation 
of solutions of As already commented in the Introduction, we shall proceed 

via F-convergence |GF75| . The reader is referred to the monographs by Attouch 
|Att84| and Dal Maso }Dal93| for some comprehensive discussion on this topic. 
Indeed, since Theorem 13.11 directly quantifies the value of the minimum to be 0, 
what is actually needed for passing to limits are F — lim inf inequalities only. We 
shall illustrate this fact by discussing the simple case of data approximation first. 

Lemma 5.1 (Data approximation). Assume (|3.5p . let £h i strongly in 
L^{0,T;Y*) being uniformly Lipschitz continuous, and i/o./i ^ J/o- Moreover, let 
Fh ■.W^'^{0,T;Y) ^ [0, oo] be defined as 



Fhiv)^ i (v^(y') + r(4-^y)-(4-A2;,y')) +x(y(0)-yo,h), 

and let Fh{yh) = 0. Then yh ~* y weakly star in W^'°°{0, T; Y) and F{y) — 0. 

Proof. Owing to Lemma 14.71 we find a (not relabeled) subsequence yh such that 
Vh ^ y weakly star in W^'°° {Q,T;Y). Hence, we have by lower semicontinuity 

(i<F{y) < liniinf|^^ {^^{y'^) + rih - Ay„) - {h,y',,)) 

+ <P{yh{T)) + <P{yo,n) ~ {Ayh{Q),yo,h) + |2/h(0) - yo,hpj 

= liminf F/i(y/i) = 0. 

Hence, F[y) = 0, y is unique, and the assertion follows from the fact that the whole 
sequence converges. □ 
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5.1. Preliminaries on functional convergence. In order to move to more gen- 
eral approximation situations, wc arc forced to discuss a suitable functional conver- 
gence notion. Wc limit ourselves in introducing the relevant definitions, referring 
to the mentioned monographs for all the necessary details. 

Recall that F is a real reflexive Banach space. Letting /„, f : Y ^ (—00,00] be 
convex, proper, and lower semicontinuous, we say that /« / in the Mosco sense 
in Y }Att841 IMos69| iff, for aU y G F , 

f{y) < hminf /„(y„) Vj/„ ^ y weakly in F, 
3y„ -> 2/ strongly in Y such that /(y) = limsup/„(2/„). 

n — 'oo 

In particular, /„ ^ / in the Mosco sense iff /„ f in the sense of F-convergence 
with respect to both the weak and the strong topology in Y. 

We will consider the situation of approximating functionals iph ■ By [Att841 Thm. 
3.18, p. 295], we have that iph V-" in the Mosco sense in Y iS tpj^ ^ -0* in the 
Mosco sense in Y*. By assuming the functionals tph to be positively 1-homogcneous, 
it turns out that the Mosco convergence tph ^ fp Y is equivalent to the Mosco 
convergence of sets — > C* in Y* which reads, 

Cn^Vn-^y* weakly in y* y*eC*, 
Wy*eC*, By; e : yl y* strongly in y*. 

Finally, we repeatedly use a lemma from |Ste06aj which we report it here, for 
the sake of completeness. 

Lemma 5.2 (Cor. 4.4, |Ste06a| ). Let p e [1,00] and fh, f ■ Y ^ (—00,00] be 
convex, proper, and lower semicontinuous such that 

f{y) < inf jliminf fhivh) ■ yh V weakly in y| V?/ G Y. 

Moreover, let yh — > y weakly in W^p{0,T;Y) (weakly star if p ~ 00), Then, we 
have that 

J{y{t)) dt< liminf / h{yh{t)) dt. 

5.2. Space approximations. We now move to the analysis of some space ap- 
proximation situation, indeed specifically tailored for the case of conformal finite 
elements. Let us enlist here our assumptions for the sake of later referencing. 



We assume to be given 

y^ C y closed subspaces such that Yh is dense in y, (5.1) 

h>0 

<t>h{y) = 4>{y) if y and (jihiy) ^ 00 otherwise. (5.2) 
i/jft : y — > (—00, 00] convex, proper, and lower semicontinuous, (5.3) 

^ph positively 1-homogeneous, (5-4) 

tph. ip in the Mosco sense in Y, (5.5) 

Hy)>^\y\'' yyeCh-Ch where = i?(V'h), (5.6) 

(■h ^ ^ pointwise strongly in y*, (5.7) 

dh uniformly Lipschitz continuous, (5.8) 

yo,/i e y^, yoji^yo- (5.9) 
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We shall mention that within the frame of conformal finite elements methods 
the subspaces Yh are obviously taken to be finite-dimensional and that the approx- 
imating functionals (f>h and iph are usually the restrictions of the functionals (j) and 
-0 on the subspace Y/j. This is exactly our choice here for 0^. In particular, one 
shall observe that ^ (/) in the Mosco sense in Y, D{d4>h) — Y^, and that 

AhV = dMy) = d(f>iy) ^ D^{y) ^ Ay yyeY^. (5.10) 



As for tph we are allowing some extra freedom (let us however remark that ([5j 
follows from (|3.5p as soon as tph is the restriction of V' to Yh since, in this case, 
Ch = CnYh). On the other hand, we are asking tph to be positively 1-homogeneous, 
namely we are considering the case of some rate- independent approximation of (jl.ip 
only. The reader is referred instead to Efendiev & Mielke [EM06j . Efendiev, 
MiELKE, Rossi, & Savare |EMRS06j . and Zanini |Zan07| for some results in the 
direction of rate-dependent approximation of rate-independent processes. 

Finally, we shall (rc)definc the approximating functionals as Fh : W^'^{Q, T; Y) 
[0, oo] as 

Fh{y) = ^ [Mv') + rhih - Auy) - (4 - y') ) + Xh{y{Q) ~ yo,h), 

where Ah = dcj)h and x/i(') — 4>h{-) + | • p. We have the following. 

Theorem 5.3 (Convergence of space approximations). Assume (|5.ip - ()5.9p and let 

Fh{yh) = 0. Then yh y weakly star in W^^°°{0, T; Y) and F{y) = 0. 

Proof. By Lemma [4.71 we find a (not relabeled) subsequence y-n ^ y weakly star 
in W^'°°{0,T;Y) and weakly pointwise. Since Fh{yh) = we readily check that 
y{t) e Yh for all t £ [0,T]. In particular, Ahyn = Ayh for all t e [0,r] owing to 
()5.10p . Hence, by lower semicontinuity, 

Q<F{y) < lin^inf(^^ (^My'h) + Vhi^h - Ayh) - {Ih.y'u)) 

+ <^(y,.(T)) + 0(yo,/.) - {Ayh{Q),yo,h) + \yh{Q) - yoA^ 
[My'h) + rh{^h-Ahyh)~{ih,y'h)) 




+ cj>h{yh{T)) ~ ^h{yhm 



= liminf Fft(y,J = 0, 



Note that the integral terms containing ip and ip* pass to the liminf by means of 
Lemma 15.21 □ 



By inspecting the proof of Theorem 15.31 (which of course generalizes Lemma 
15. ip . one realizes that, whenever the weak-star precompactness in M^^'°°(0, T; Y) of 
the sequence yh is assumed, the convergence statement holds more generally in the 
case Fh{yh) 0. Namely, by directly asking for the above-mentioned compactness, 
one could consider the convergence of some approximated solutions yh such that, 
possibly, Fh{yh) > 0. We rephrase this fact in the following statement. 
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Lemma 5.4 (r-liminf inequality for Fh). Assume ((5?T|) -([53 )) . ((53)) - ((ST)) . and 
dSH). Then, 

F{u) < inf i limM Fhiyh) ■ Vh -> V weakly star in VFi'°°(0, T; Y) 

Note that the homogeneity of iph ■, the uniform convexity of 0^ , and the Lipschitz 
continuity of play no role here. 

Finally, again by looking carefully to the proof of Theorem 15 . 31 one could wonder 
if the requirement on the Mosco convergence of could be weakened. Indeed, 
what we are actually using is only that 

?/> < r - lim inf and V* < T - lim inf (5.11) 

with respect to the weak topologies of Y and Y* . respectively. On the other hand, 
in our specific situation, [Ste06a[ Lemma 4.1] entails that (|5.11|) and the fact that 
■0/1 ^ t/i in the Mosco sense in Y are equivalent. 

This observation motivates once again the belief that Mosco convergence is the 
right frame in order to pass to limits within rate- independent problems. For the 
sake of completeness, let us recall that a first result in the direction of the approx- 
imation of the play operator {Y Hilbert and A coercive on F) under the Hausdorff 
convergence of the characteristic sets C,* = D{il'fJ is contained in [Kre961 Thm. 
3.12, p. 34] whereas the extension of this result to the more general situation of 
Mosco converging sets as well as some application to parabolic PDEs with hystere- 
sis is discussed in jSte06b| . More recently, Mielke, Roubi'cek, <k Stefanelli 
[MRS06| addressed in full generality the issue of F-convergence and relaxation for 
the energetic solutions of rate-independent processes. An alternative convergence 
proof in the specific case of convex energies is obtained by means of the Brezis- 
Ekeland approach in [Ste06a| . 



6. Time-discretization 

Assume now we are given the partitions P„ = {0 = t° < < • • • < t^" = T} 
and denote by r,* = — t^^^ the i-ih time-step and by t„ = maxi<i<jv„ the 
diameter of the n-th partition. No constraints are imposed on the possible choice 
of the time-steps throughout this analysis besides t„ — *■ as n — > oo. Moreover, let 
the parameter 9 S [1/2, 1] be given. 

In the following we will make an extensive use of the following notation: letting 
V = {v^, . . . , w^") be a vector, we will denote by t;„ and TJ„ two functions of the 
time-interval [0, T] which interpolate the values of the vector v piecewise linearly 
and backward constantly on the partition P„, respectively. Namely 

z?„(o) = v„it) = ^i^ity + (1 - 

W„(0) = i;", v^{t)^v\ for i e i = l,...,iV„ 

where 

7; W - - for t e z = 1, . . . , iV„. 

Moreover, we let Sv"^ ~ (w* — v^^^)/Tn for i = 1, . . . , iV,i (so that v'^ — Svn) and 
denote by ve the vector with components I'l = Ov"^ + (1 — d)v^^^. 
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Recall that t e W^^°"{Q,T;Y*) and yo e 5(0). We shall be concerned with the 
so-called 0-scheme for problem (|l.ip : 

( ^n ~yn' \^ ^ ^(gy. ^ _ 3 ^(^^j^ ^ _ 

for z = l,...,7V„, (6.1) 

vl - yo. (6.2) 

One usually refers to the latter as backward or implicit Euler scheme for the choice 
9 = \ and as Crank-Nicholson scheme for 6* = 1/2. 

Owing to the above-introduced notation, the latter scheme can be equivalently 
rewritten as 

5^(2/;-2/r')+^<e9^(iL(^) for i = l,...,iV„, 2/°=yo. (6.3) 



Clearly, the 6'-scheme ()6.3p is rate-independent. Namely, no time-step appears in 
(|6.3|) and the choice of the partition affects the solution via the values of the load £ 
only. In this concern, our focus on variable time-steps partition could be simplified 
by considering proper rescaled loads £ instead. We shall however keep up with it, 
especially in order to underline the possibility of adapting the partition according 
to some a posteriori analysis (see Subsection 16. 7p . 

Before moving on, let us comment that, for all n, the latter scheme as a unique 
solution. Indeed, given y\^^ £ C, it sufhces to (uniquely) solve iteratively the 
incremental problem 

e AYgmmy^y[0^{y) - {e{tl,) - (1 - 0)Ayl-\y) + i,{y - y^-')). (6.4) 

Note that, since y^^^ € C, the functional under minimization turns out to be 
uniformly convex. Hence, by (j3.4[) . the minimum problem has a unique solution. 
In particular, exactly as in Lemma [Ol we have the following. 

Lemma 6.1. yl^ £ C for all i = 0,1, ... , Nn. 

A crucial observation is that, as in the continuous case, the discrete trajectories 
show some sort of stability as well. 

Lemma 6.2 (Stability of the discrete trajectories). We have that 

y; e Argmin^gy.(00(y) - (£(ij,,) - (1 - 0)Ayl-\y) + ij{y - yjj) 

for i = l,...,Nn. (6.5) 
In particular, if 9 = \, then j/^ € S{t\^). 



Proof. From the incremental formulation (|6.4p and the triangle inequality for -0 we 
get that, for all y G F, 

e^t^iVn) - (iitle) - (1 - 0)Ayl^-\yl^) + ^{yl^ - f-') 

< 9<P{y) - (iitl,) - (1 - 9)Af-\y) + ^{y - y^') 

< 9,p{y) - {£{tl,) - (1 - 9)Ayl-\y) + 0(y - ylJ + - yl~'), 

whence the assertion follows. □ 

Again as in the continuous case, we readily check that 

relation (113]) holds iff £{tl g)-Ayl^j)eC*. (6.6) 
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6.1. The discrete variational principle. We shall now present a discrete version 
of the variational principle of Theorem 13.11 

We define L^'^iy, z) : Y x Y ^ [0, oo] as 

Li^\y, z)^yj (^^) + r {e{tl,) - A{0y + (1 - e)z)) 

eiti,)-Aiey + ii-0)z),y—^ 

' n 

and the functionals : [0,oo] as 

F^{yl ...,y^-) = Y^ <i4^{y^, Vn') + x{yl - yo). 
1=1 

Lemma 6.3 (Discrete variational principle). . . . , y^") solves 
F«(yO,...,y^)=0 = mini^«. 

Proof. Analogously to the continuous case, we have that, for all i = 1, . . . , A^„, 

di^ (<52/;) + Ayl, 9 iff y^') - 0, 

and yl^ ^ yo iff x(y" - yo) = 0. □ 

Let us observe that the functional Ff^ is convex and lower-semicoiitinuous. More- 
over, by the homogeneity of ip (see p.2p ). F.^ is actually independent of the time 
steps. In fact, we have 

Jv„ / 



- (eitlg) - Ayl„ yl, - yr 1) j + xivl - 2/o). 

The idea of dealing with time-discretizations via a discrete variational principle 
closely relates our analysis to the theory of so-called variationals integrators. The 
latter are numerical schemes stemming from the approximation of the action func- 
tional in Lagrangian Mechanics. By referring the reader to the monograph |HLW06| 
and the survey |MW01| , we shall restrain here from giving a detailed presentation of 
the subject and limit ourselves to some (necessarily sketchy) considerations. Let- 
ting {t,y,p) £ [0,T] X R™ X M™ C{t,y,p) denote the Lagrangian of a (finite 
dimensional, for simplicity) system, the Hamilton principle asserts that the actual 
trajectory 1 !■ y{t) of the system minimizes the action functional 

y^ f C{t,y{t),y'{t)) dt 



Jo 

among all curves with prescribed endpoints, thus solving the Lagrange equations 
dy-C ~ -^dp.C = for i = 1, . . . , TO. (6.7) 

Hence, a natural idea is that of deriving numerical schemes for Lagrangian mechan- 
ics by applying some quadrature procedure to the action functional, i.e. discretizing 
Hamilton's principle. The resulting discrete schemes show comparable performance 
with respect to other methods but generally enjoy some interesting extra (and of- 
ten crucial) properties such as the conservation of suitable quantities |LMOW04] . 
Variational integrators have been intensively applied in finite-dimensional contexts 
and, more recently, to the situation of nonlinear wave equations |MPS98j and non- 
equilibrium elasticity [LMOW03| . 
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The present analysis may bear some resemblance to the above-mentioned theory. 
Indeed, the formulation of the 6'-scheme in the case 9=1/2 stems exactly from the 
midpoint quadrature of the functional F as 

= ^{y^^y^-^)+r{i{-t^nS,2)~A{y^ 

where y is taken to be piecewise affine on the partition P„ . 

On the other hand, our focus here is quite different. First of all, we are not 
dealing with the Hamilton principle (cndpoints are not fixed) as we are not aimed 
at solving the Euler-Lagrange equations for F (i.e., solve (|6.7p ). Secondly, we are 
specifically interested at infinite-dimensional situations, namely PDEs. Finally, the 
only choice of 9 which is directly related with a quadrature of F is = 1/2 and we 
are not considering higher order schemes. 

Before closing this discussion, let us mention that some F-convergence techniques 
have been recently exploited in the (finite-dimensional) frame of variational inte- 
grators by MuLLER & Ortiz |MO04j (see also |MM04j ). In this same spirit, we 
are here providing F-convergencc results in infinite dimensions instead. 

6.2. Stability of the 6l-scheme. It is known since Han & Reddy |HR95[[HR99] 

that the choice 9 < 1/2 in (|6.3p leads to an unconditionally unstable scheme and 
that, on the contrary, for 9 G [1/2, 1] the 0-scheine is stable in i/^(0,T;y) when Y 
is a Hilbert space and the partitions arc chosen to be uniform. 

Here we shall provide an alternative stability proof by taking into account the 
Banach frame. 

Lemma 6.4 (Stability). Assume p.5p and let 9 G [1/2,1]. Then, the solution to 
the 9-scheme ()6.3p fulfills 

\\y'n,e\\L-iO,T:Y) < W iO,T;Y') = 1 OT 9 =\. (6.8) 

Moreover, for constant time-steps, 

\\y'n,0h-(O,T;Y) < J„ H lU^ (o,r;y ♦ ) tf ^ < 9 < 1. (6.9) 

Our argument coincides with that of |MT041 Thm. 4.4] in the case of Euler, 
i.e. 9 = 1 and it is an extension of the latter for the case 1/2 < < 1. Here, 
we do not play with the variational inequality by choosing suitable tests but use 
the scalar relations L^'*(2/^, y^~^) = instead (this makes however no substantial 
difference since the latter scalar relations are exactly the outcome of the test on the 
variational inequahty in jMT041 Thm. 4.4]). 

The stability proof for the Crank-Nicholson scheme 9 = 1/2 is quite differ- 
ent from former arguments and stems as a direct outcome of our variational ap- 
proach. Let us mention that, unlike the classical parabolic situation, here the 
Crank-Nicholson scheme is indeed unconditionally stable. The reason for this fact 
is the rate-independence of the problem, namely the degenerate character of the 
evolution (no relaxation time). In both cases 9 = 1 and 9 = 1/2, the stability 
constant 1/a is sharp (see Lemma FZ)) . 
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We complement this analysis by providing the stability for the 0-scheme for 
1/2 < < 1 in the case of constant time-steps (likely with a non-optimal, although 
explicit, stability constant). 

Proof. Let us prove the stability of the Crank-Nicholson scheme 9 — 1/2 first. For 
this aim, it suffices to recall that 

0==Fy^(yJl,...,y„^") 

= (iiffj + r{l. - Ay„) - (?„ - Ay„, y;) ) + x(y„(0) - yo). 

Hence y„ minimizes the functional F where £ is replaced by The stability 
estimate follows from Lemma 14.71 

Let us now move to the case 1/2 < 6' < 1. Relation (|6.5p applied at level i — 1 
for some i = 2, . . . , Nn along with the choice y — yl^ entails that 

< 00(2/;) - (nC,i) - (1 - e)Ayl~',yC) + i^{yl, - ^r^), 

where the extra-term 6(j3[yl^ — y\^^) is obtained from the fact that (j) is quadratic. 
Hence, we have that 

e4>{yn-yn') + 0c^{y'n')-04>{Vn) 

< {^KD^yn' - yn) + (A(yr' - yr'), - Vn') + i'iVn - Vn') 

+ {l-0){Af-\yl-yl-') 

= {^itn,e)^yn' - y'r) + (1 " ^) (^(y^' - yr'), y; - Vn') + i^ivl ~ Vn') 

-{1 - 0)(0(yr') + m - yr') - ^(y;)) , 

so that 

Hen) + 4>ifn-') - 0(y^) < - (eKj), el,) + {e-i) (Aejr\ + V'(e^), (6.10) 

where we have used ejj ~ yl, — y]^^ in order to shorten notations. 
Next, from Lj, ^{yl^, yl^^) = for i = 1, . . . , Nn, we obtain that 

= V'«)-(^(t;,J-A<e,<) 

= ^(ej,) - {i{tle),el) + 0{Hyn) + HO - Hy'n')) 

- (i-0)(0(yr') + 0(ej.)-0(y;)). 

In particular, we have checked that 

^{el,) + Hyn) - 0(yr') + (20 - men) = (^(t^e), ■ (6.11) 
We take the sum between the latter and (|6.10p and get that 

2eHel,) < {litis) - £(Ce + (0 - 1) {Ael-\el) , 

or, equivalcntly, 

{Aelg,el) < (e{tl,) - i{tl'g'),el,) . (6.12) 
Now, if 6* = 1. we conclude that 

K\<\<l^mio)-eitn:e)\*^ 

and the assertion follows. 
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In case 1/2 < 9 < 1 and for a constant time-step partition, one proceeds from 
(|6.12p by computing 

'2 



-r™ II CII L°° (o,T;y* ) Y - V^RJ 

>(?(A4,e:,) + (l-0)(Ae:r\e:,) 

= [29 - 1) {A^, 4) + {l-9) (A(e; + ejr'), ej,) 

> 2(20 - l)</>(e;) + (1 - 0)(^(e;) - 0(6^^)) 

Note that the coefficient [29 — 1) is strictly positive as 6* > 1/2. By using the fact 
that ?/° = yo £ '5(0) (recall (|4.5p ). we readily check that 

^{el) + cjy{yl) ~ (£(0), yl) < (b{yl) - (£(0), yl) + ^(e^ ) 

and, by adding the latter to (|6.1ip for i = 1 we have 

20</>(e,\) < {e{tlo) - m,4) < rn\K\\L^io,T;Y')^V^)- (6.13) 
Let us define 



2(20-1) _ 1 [2^ Ci 

so that, owing to (|6?T3l) . ((6TT3ll . and using the fact that 2(26* - 1) < 20, 
(Co + l)a^-ati < Cia, for i = 2,...,iV„, 

^ T^V?"^""'^""'"'^^^*^ 2(20-1) 
C'l 

— TT - <^2- 

Now, since (Cq + 1)C| — C1C2 — C|, we easily prove by induction that Ui < C2 
and the assertion follows. □ 



6.3. Convergence. We shall prove the weak-star W^'^{0,T;Y*) convergence for 
the 0-method. This result has to be compared with that of Han & Reddy [HROOi 
Thm. 3.4] where the uniform convergence of the backward constant interpolations 
is obtained. Our result is weaker than that of [HROOi Thm. 3.4] since we are not 
providing strong convergence. On the other hand, we believe our half-page proof to 
be possibly more transparent than the long argument developed in |HROO| . Let us 
moreover mention that in the Hilbertian case and for A coercive on Y, the strong 
convergence in W^'P{0,T;Y) for all p < 00 of the Euler method = 1 has been 
proved in |Kre96[ Prop. 3.9, p. 33]. 

Theorem 6.5 (Convergence for the 0-method). Assume (|3.5p and let 

F^{y°,...,y^") = 0. Then, y„ ^ y weakly star in W^°°{0,T-Y) where 
F{y) = = minF. 
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Proof. Owing to Lemma 16.41 we can extract a (not relabeled) subsequence such 
that Vn ^ y weakly star in W^ "^{0, T; Y), hence weakly pointwise in Y. Moreover, 
we clearly have that both y„ and y„ g converge at the same limit weakly star in 
i°°(0, T; Y). Finally, we directly check that Z„,e £ strongly in L°°{0, T;Y*). By 
observing that, since 6 > 1/2, 

< {AWn + (1 - d)yn-'),Syl) = m.) + (20 - " yn') " 'PiyV) 

> <^(y:j-0(yri). 



we compute that 



= F:(y::,...,y-") 

[i^iyn) +'>P*(^n,e - Ay^ g) - {£n,e,yn)) 







> 

+ (t>{yn{T)) - + x(yn(0) - yo) 

^ 

4>{yn{T)) + </)(yo) - (Ay„(0), yo) + |yn(0) - yo 



2 



Finally, it suffices to pass to the liminf above as n oo and exploit lower semi- 
continuity and the stated convergences in order to obtain that F{y) < 0. Hence, 
by Theorem 13.11 and Corollary 14. 6( y is the only solution to (|l.ip and the whole 
sequence y„ converges. □ 

We shall mention that the separability assumption for Y is not crucial and could 
be weakened. Indeed, in case Y is not separable, one simply has to pass to limits 
as Tjn y weakly in W^'P{0, T; Y) for some p S [1, oo). 



6.4. The functional controls the uniform distance. Wc shall reproduce at the 
discrete level the results of Subsection 14.61 We begin by showing how to possibly 
control the uniform distance of two vectors by means of the discrete functional F^. 

Lemma 6.6 (Uniform distance control via F^). Let the vectors u = (m°, . . . , u^"), 
V = {v^ , . . . , v^") G y^n+i 6e given. Then, 

77(1 — 77) max (j)(u^ — v^) 

l<i<Ar„ 

< r,F^{u', . . . ,^.^") + (1 - v)F^Av', . . . V77 e [0, 1]. 

Proof. This proof follows the same lines of that of Corollary 14.41 Let 1 < i < Nn 
be fixed and define F^'^ : [0, 00] as 

F^-\y',...,f) 

i 

= Y.<Ln'iy',y'-') + x{y°-yo) 

i 

= 5]V(y'' -2/^-') + V'*(^(D - Ay^) - {m,e)^y' -y'^') 

i 

+ <j>{y') + {29 - 1) ^ 0(2/^- - + <j>{y„) - (V, yo) + |y" - yoP- 
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Then, clearly y = . . . , y*) i— > G^'^{y) = F^'^{y) — (j){y^) is convex. Hence, letting 
w = r]u + {1 — r])v for some 77 G [0, 1] we have that 

< F^'^w) 

< v{G'/{u) + + (1 - v){G':fiv) + - - v)'P{u^ - v'), 

whence the assertion follows. □ 

Again, note that the latter lemma controls the uniform norm of the distance only 
if the stronger (j3.5p is required. The following corollary of Lemma [6.61 will be the 
starting point for some possible a posteriori error control procedure (see Subsection 
Ml- 

Corollary 6.7 (Uniform distance from the minimizer). Let F^^{y'^, . . . ,y^") = 0. 
Then 

^max ^{y'~v')<F^{v°,...,v^") V(v°, . . . , t;^") € 

Moreover, wc re-obtain a proof of the uniqueness of the solution of the 0-method. 

Corollary 6.8 (Uniqueness of the minimizer). Assume p.Sp . Then, there exists 
at most one y = {y^ , . . . , y^") such that F^^{y) = 0. 

6.5. The generalized 6'-method. Although minimizers of F^ and solutions of 
the 6'-scheme (|6.3p coincide, minimizing sequences of F,j need not solve (|6.3p . This 
extra freedom allows the minimization formulation to capture the convergence of 
some generalized 0-method, where the relations in ()6.3p are not solved exactly but 
rather approximated. Namely, we shall look for vectors m„ = [u^, . . . such 
that 

Fn{un) — * as n — ^ 00 
instead of F^{un) = for all n e N. 

From the computational viewpoint, note that the ^-scheme consists in solving 
Nn nonlinear equations in one unknown each while checking for stationarity for 
F^ implies the solution of a tridiagonal system of Nn + 1 nonlinear equations 
with (up to) three unknowns each. This entails in particular that minimizing F^ 
instead of solving (|6.3p could be of a scarce interest if one is merely concerned in 
reproducing the ^-scheme with no error. On the other hand, the issue of solving 
up to some tolerance turns out to be particularly relevant whenever one is aimed 
at implementing an optimization procedure for the solution of (|6.3p . Indeed, one 
should be prepared to run the algorithm (some descent method, say) until some 
given tolerance is reached. 

Our starting point for a possible convergence analysis of the generalized 0-method 
is the following classical error control result. 

Theorem 6.9 (Mielke & Theil |MT04j ). Assume (P^]) . Then, y„ ^ y uniformly 
and F(y) ~ 0. In particular, 

max |(y„ - y)it)\ <C,^ (6.14) 

tG[0,T] 

where Ce depends only on data and is independent of n. 

More precisely, in |MT04j solely the case of the Euler scheme 6* = 1 is discussed. 
However, an easy adaptation of the argument entails the result for 6 G [1/2, 1) as 
well. 

By explicitly comparing the minimizing sequence u„ = (m°, . . . ,m^") with the 
corresponding solution (?/°, . . . , ?/^") of the 0- method, we have the following. 
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Theorem 6.10 (Convergence for the generalized ^-method). Assume (|3.5p and let 

F^{u'^, . . . ,u^") 0. Then, u„ — *■ y uniformly, where F{y) = 0. In particular, 

max |(^2„-2/)(t)|<CeV^+ . (6.15) 



*G[o,T] \a. 

Proof. We have that 

max |(y-u„)(i)| < max |(y - y„)(t)| + max |(y„ - u„)(i)| 

< CeVT^ + rnax - y'J 

1<1< Nn 



< Cey^r + ( - max 0(uj, - y'J 
a i<i<iv„ 



1/2 



and we conclude by applying Corollarv l6.7l □ 



6.6. A posteriori error control. Let us now exploit both CoroUarv 14. 51 and The- 
orem 16.101 in order to provide some possible a posteriori estimates of the approxi- 
mation error by means of solutions u„ of the generalized ^-method described above. 

Lemma 6.11 (A posteriori error control via F^). Assume ()3.5p and let 

F^{u^, . . . , w^") ^ for some s > and F{y) = 0. Then, 

max \(un — y)it)\ ^ tII where 2r = max{l, s}. 
te[o.T] 

Lemma 6.12 (A posteriori error control via F). Assume (j3.5p and let F(u„) ~ 
for some s > and F{y) = 0. Then, max4g[o,T] K^n ~ y){t)\ ~ Tn^- 

We are also in the position of proving the weak-star convergence of the time 
derivatives of solutions u„ of the generalized ^-method by comparing them with 
the corresponding derivatives of the exact solution of the 0-method. 



Lemma 6.13 (Improved convergence for the generalized 6'-method). Assume ()3.5p 
and let . . . , it^") ~ r^. Then, u„ is equibounded in W^^°°{0,T;Y). In 

particular, u„ — > y weakly star in W^'°°{0,T;Y). 

Proof. Let (y°, . . . , y^") be the solution of the 6'-scheme. By exploiting Lemma 
16. 7( we check that 

\<-<'\ < \<-yn\ + \yn~yn~^\ + \yn~^ -<~'\ 

,1/2 

< 



r;l|y;i|L-(0,T;r) + 2 (^^F„««, . . . , «^) 



The uniform bound on ||un||vi/i'°°(o,T;F) follows by dividing the latter by r^, taking 
the maximum as 1 < z < Nn, and recalling Lemma [6.41 □ 



6.7. Adaptivity. Assuming (|3.5p . the above introduced a posteriori error estima- 
tors can be exploited in order to develop an adaptive strategy. In particular, the 
error control in the uniform norm up to a given tolerance tol > 

max |(y - y„)(t)| < tol 
*e[o,T] 

for some piecewise approximation yn with x(yn(0) — yo) ^ « tol^/4 can be inferred, 
for instance, by choosing time steps in such a way that 

*" , / ^ / s w / N N Oi tol^ 
L(i,y„(t),y;(i))< — — . 

Namely, by uniformly distributing the error along the partition. 
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Alternatively, one could develop an adaptive strategy by considering just com- 
puted quantities at the discrete level by asking for 

i Td i/ i ^ atol^ r ^ tol^ 

<i„ ) < T^TTTT for r„< 



327V„ - 16C3 



and exploiting Theorem 16. 101 



7. SPACE-TIME APPROXIMATIONS 

Let us combine the results of the previous sections (and use the corresponding 
notation) in order to state and prove a result on the convergence of full space-time 
approximations. Our results have to be compared with the former convergence 
analysis by Han & Reddy [HR99j . Our approach leads to a convergence proof 
with respect to a weaker topology. However, it is on the one hand slightly more 
general (some assumptions on the spaces and the functionals, sec (H1)-(H2) |HR99[ 
p. 264], are not required) and on the other hand has a much simpler proof. 

Theorem 7.1 (Convergence of space-time approximations). Assume ()5.ip - (|5.9p . 

9 € [1/2, 1], define i^;^^' ^) - Y xY ^ [0,oo] as 

<Uy, ^) - i'h (^) + iUtle) - A„{ey + (1 - 0)z)) 
- (e„{tlo) - A,,{0y + (1 - 0)z), y—^^ , 
where Ah = d4>h, and let the functionals : y^^+i [0, oo] he defined as 

1=1 

where Xh{-) = (t>h{-) + I ■ P (note that D{Fl^) C Y^"+^). Finally, let 
F„j,(yO,...,yf") = 0. We have: 

(a) yn,h —* yh weakly star in W^'°°{0, T; Y) as (n, h) (oo, h) and Fh{yh) = 0. 

(b) yn,h y-n weakly star in W^'°°{Q,T;Y) as {n,h) {n,0) and F^{y„) = 0. 

(c) yn^y weakly star in W^^°°{0,T;Y) as (n, 0) {oo,0) and F{y) = 0. 

(d) yu y weakly starW^^'^{Q,T]Y) as (oo,/i) (cx),0) and F{y) ^ 0. 

(e) ynji —* y weakly star in W^'°°{0, T; Y) as (n, h) —f (oo, 0) and F{y) = 0. 

The thesis of the Theorem is illustrated in Figure 1. In particular, we aim at 
showing that the space (or data) and time-limit can be taken in any order. Note 
that Limit (c) has been already checked in Theorem 16.51 and that the very same 
argument yields Limit (a) as well (recall that Y^ is closed). Moreover, Limit (d) is 
discussed in Theorem 15. 31 So what we are actually left to check arc Limits (b) and 
(e) only. 

Proof. Limit (b). The assertion follows once we check that, for alH = 1, . . . , ^n, if 
VrTh ^ yr7^ weakly in Y, one has the weak convergence ynh^Vn ^ ^^U- Recall 
that 

ylh e Argmmy^Y{sMy) - (hit^e) " (1 " 0)Ahyl-h^y) + My - Vnl)) 

= Argmin^^y,^ (^^(y) - (4(iLe) - (1 - 0)Ay^;,,y) + My - yl~l)) ■ 
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d 




• 



(c>D,0) n 
Figure 1: Convergences for space-time approximations 



Hence, since we have (|5.6p . the sequence is weakly precompact and, up to the 
extraction of a (not relabeled) subsequence, Un h ^ V weakly in Y . Let us prove 
that y solves the incremental problem (|6.4p . Indeed, we have that 

< i^'Xy,2/r') 

< liminf (MVn,h ' vlZl) + rd^hi^nfi) - M^n.Ke) 

= lim inf L^';^ (y; , y'-,] ) = 

where we have used the Mosco convergence in (jS.Sp and the pointwise convergence 
of ih (|5.7p . Since the only solution of (|6.4p is y^, we have that y = and the 
whole sequence converges. 

Let us mention that, if the functionals Tph are uniformly linearly bounded (which 
is quite common in practice), then one could prove the latter convergence to be 
actually strong: namely y^^jl — > yl^^ strongly in Y implies the strong convergence 
Vnh Vn- Indeed, let Wh and Wh be such that Wh — y^ — > strongly in Y, 
■>Phiwh - yl, f^) ^ 0, w/i e Yh and Wh ~- Wh ^ strongly in Y. Then 

eHylh) - {e{tle) - (i - o)A^fJ,yi,,) + MvIh ~ vTm) 

< e4>h{wh) - {iitn,g) - (1 - 9)Ayi-^,Wh) + Mwh - Wh) + Mwh - 2/L J. 

If iph are uniformly linearly bounded above, then iphiwh — Wh) — > with /i ^ 0. 
Whence, by passing to the lim sup in the latter, we check that limsup,j_o (piVn h) — 
4>iyn), which, together with lower semicontinuity gives (^(y^ 4>{yn) and the 

strong convergence follows from the reflexivity of Y . 

Limit (e). Lemma I6.4| the uniform Lipschitz continuity of (|5.8|) . and the 
initial datum convergence (|5.9|) . entail that y„_/j are uniformly Lipschitz continuous 
as well. Hence, by extracting a (not relabeled) subsequence, yn,h y weakly star 
in W^'°°{0,T;Y). In order to check that y solves (jl.ip . let us remark that, being 

ln,hfi^l Strongly in L^{Q,T;Y*) 
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and that, by [Ste06a[ Cor. 4.4] 







T 

i;*{£ - Ay) < liminf / ^K^nMfi - ^y„,h,e), 



and compute that 
< F{y) 




+ 0(yri,h.(r)) + (/)(yojO - (Ay„_/,(0), yo,/i> + IVnA^) - Va^l j 

and we have that F(y) = 0. □ 

We shall conclude by briefly mention some further results which can be obtained 
by suitably adapting to the current fuUy-discretized situation the arguments devel- 
oped above for timc-discretizations. Firstly, in the same spirit of Lemma 16.61 one 
could consider the possibility of estimating the distance of a vector from the mini- 
mizcr of ^ by means of the functional itself. Secondly, the use of Corollary 14.51 
would entail the possibility of an a posteriori error control an some adaptive strat- 
egy along the lines of Subsection 16 . 71 could be considered. Finally, by relying on the 
known convergence estimates for full space-time discrctized problems |HR99j one 
could obtain a convergence and an a posteriori error control result for some gener- 
alized space-time approximated problem where F^^ are not exactly minimized and 
one considers minimizing sequences instead (see Subsection 16. 5p . We shall develop 
these considerations elsewhere. 
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